Permeative flows in cholesteric liquid crystals 
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We use lattice Boltzmann simulations to solve the Beris-Edwards equations of motion for a 
cholesteric liquid crystal subjected to Poiseuille flow along the direction of the helical axis (per- 
meative flow). The results allow us to clarify and extend the approximate analytic treatments 
currently available. We find that if the cholesteric helix is pinned at the boundaries there is an 
enormous viscosity increase. If, instead, the helix is free the velocity profile is fiattened but the 
viscosity is essentially unchanged. We highlight the importance of secondary flows and, for higher 
flow velocities, we identify a flow-induced double twist structure in the director field - reminiscent 
of the texture characteristic of blue phases. 
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Liquid crystals are fluids, typically comprising long 
thin molecules, where subtle energy - entropy balances 
can cause the molecules to align to form a variety of 
ordered states [1,2]. In nematic liquid crystals the 
molecules tend to align parallel giving a state with long- 
range orientational order. This is usefully described by 
the director field ft, the coarse-gained, average, molecu- 
lar orientation. In a cholesteric or chiral nematic liquid 
crystal n has a natural twist deformation in the direction 
perpendicular to the molecules. Examples of cholesteric 
liquid crystals are DNA molecules in solution, colloidal 
suspensions of bacteriophages [3], and solutions of ne- 
matic mixtures such as E7 with chiral dopants which are 
widely used in display devices. 

Liquid crystals exhibit both an elastic and a viscous 
response to an external stress. Coupling between the 
director and the velocity fields - known as back-flow ~ 
leads to strongly non-Newtonian flow behaviour. A par- 
ticularly striking example in cholesterics is permeation. 
When a cholesteric liquid crystal is subjected to an im- 
posed flow in the direction of the helix axis, its viscosity 
can increase enormously (by a factor ~ 10^) when the 
isotropic to cholesteric transition is reached [4-6]. 

An explanation of permeation was given by Helfrich 
[4]. If the director orientation is fixed in space, due for 
example to anchoring effects at the wall, any flow along 
the helix must be linked with a rotation of the molecules. 
This leads to an energy dissipation far larger than that 
due to the usual molecular friction and hence a much 
enhanced viscosity. 

Balancing the dissipation from the director rotation 
with the energy gained from the pressure gradient along 
the capillary, Helfrich argued that the usual parabolic 
velocity profile is replaced by plug-like flow, with a con- 
stant velocity across the capillary [4]. To satisfy no-slip 
boundary conditions, however, the velocity must fall to 
zero at the edges of the sample. De Gennes and Prost 
[1] argue that this occurs over a length scale ~ p where 
p is the pitch of the helix. Rey has extended these ideas 
to linear and oscillatory shear [7]. These calculations all 
assume small forcing, so that the director field is not de- 
formed by the helix. There is an interesting suggestion 
in an early paper by Prost et al. [8] that an increase in 



forcing might stabilize a modulation in the director field 
in a direction perpendicular to the flow. 

There are few quantitative experiments on cholesteric 
rheology, mainly as it is difficult to obtain single domain 
textures or a uniformly pinned helix [5,6,9]. Porter et al. 
showed firm evidence for a high cholesteric viscosity at 
low forcing [5]. This dropped to values close to those of 
a normal liquid crystal as the shear rate was increased. 

Helfrich's explanation of permeation is widely accepted 
but many questions remain. There is confusion in the lit- 
erature about whether the director field must be pinned 
in some way to obtain a permeative flow. It is interesting 
to ask whether distortions in the director field, induced 
by the flow, alter the permeation. Does permeation per- 
sist beyond the regime of low forcing and what replaces 
it for larger forcing ? Finally, how does the interplay be- 
tween the width of the boundary layer and that of the 
channel affect the flow? 

Given the importance of cholesteric liquid crystals in 
optical devices and biological DNA solutions and the 
ubiquity of permeative flow in the theory of layered liquid 
crystals such as cholesterics and smectics [10] it is use- 
ful to develop a robust numerical method for solving the 
equations of cholesteric hydrodynamics, thus allowing us 
to probe questions that cannot be answered analytically. 
Moreover the advent of micro-channel technology means 
that quantitive experiments are likely to become feasible. 

Therefore in this paper we solve numerically the hydro- 
dynamic equations of cholesterics and hence clarify some 
of the outstanding questions about permeation. A major 
conclusion is the importance of boundary conditions. If 
the helix is pinned at the boundaries, e.g. by wall ir- 
regularities, the apparent viscosity increases by orders of 
magnitude. If the boundaries are free, the velocity profile 
is flat but the helical structure is free to drift along the 
flow direction and the viscosity is much smaller. As the 
velocity of the system is increased shear forces induce 
a signiflcant deformation of the initial helix. First the 
cholesteric layers are bent into chevrons. Then at higher 
forcing a doubly twisted texture is formed, the initial de- 
formation being accompanied by a flow-induced twist in 
the perpendicular direction. 

We consider the formulation of liquid crystal hydro- 
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dynamics given by Boris and Edwards [11]. The equa- 
tions of motion are written in terms of a tensor order pa- 
rameter Q which is related to the direction of individual 
molecules, h, by Qa/J = {nafi0 — 3<^a/3) where the angular 
brackets denote a coarse-grained average and the Greek 
indices label the Cartesian components of Q. The tensor 
Q is traceless and symmetric. Its largest eigenvalue, ^q, 
< q < 1, describes the magnitude of the order. 

The equiUbrium properties of the Uquid crystal are de- 
scribed by a Landau-de Gennes free energy density. This 
comprises a bulk term (summation over repeated indices 
is implied hereafter). 
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and a distortion term, which for cholesterics is [1] 
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where K is an elastic constant. The tensor e^^^ is the 
Levi-Civita antisymmetric third-rank tensor, Aq is a con- 
stant and 7 controls the magnitude of order. The an- 
choring of the director field on the boundary surfaces 
is ensured by adding a surface term proportional to 
{Qap — Qq/j)^: with chosen in such a way that the 
director has the desired orientation at the boundaries. 
The equation of motion for Q is [11] 



(a* +u- V)Q-S(W,Q) = rH 



(3) 



where F is a collective rotational diffusion constant. The 

first term on the left-hand side of Eq. (3) is the mate- 
rial derivative describing the usual time dependence of 
a quantity advectcd by a fluid with velocity u. This is 
generalized for rod-like molecules by a second term 

S(W, Q) = (CD + n){Q + 1/3) + (Q + I/3)(CD - fi) (4) 
- 2C(Q + I/3)Tr(QW) 

where Tr denotes the tensorial trace, while D = (W + 
W"' )/2 and ft = (W — W^)/2 are the symmetric part 
and the anti-symmetric part respectively of the velocity 
gradient tensor Wap = dpUa- The constant ^ depends on 
the molecular details of a given liquid crystal. The term 
on the right-hand side of Eq. (3) describes the relaxation 
of the order parameter towards the minimum of the free 
energy. The molecular field H is given by 
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The three-dimensional fluid velocity, u, obeys the con- 
tinuity equation and the Navier-Stokes equation, 

p{dt + Updfj)Ua = dij{l\.al3) + Vd)3{daUf3 + d^Ua (6) 
-I- (1 - ^dpPo)d^U^5al3), 



where p is the fluid density and rj is an isotropic viscosity. 
The stress tensor 11^0 necessary to describe liquid crystal 
hydrodynamics is: 



Ha/? = 
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Pq is a constant in the simulations reported here. 

The differential equations (3) and (6) are coupled. Un- 
less the flow field is zero (u = 0) the dynamics given by 
Eq. (3) are not purely relaxational. Conversely, the or- 
der parameter field affects the dynamics of the flow field 
through the stress tensor. This is the back-flow coupling. 
To solve these equations we use a lattice Boltzmann al- 
gorithm. Details and validation of using this method to 
solve the Beris Edwards model were given in Ref. [12,13]. 

We consider a cholesteric liquid crystal which is sand- 
wiched between two plates a distance L apart along the 
z-axis. The axis of the cholesteric helix lies in a direc- 
tion parallel to the plates which we shall take as the y- 
axis. The primary flow is also along y. and is imposed 
via a pressure gradient. In the steady state, however, 
there can be secondary flow so that the modelling must 
be fully three-dimensional. This geometry is the one for 
which the permeation mode is expected. The elastic con- 
stants and viscosities are taken to have typical values for 
a cholesteric and wc consider channel widths L ^ fim. 
Typically, a cholesteric pitch was discretised by 64 lat- 
tice points, while 10^ iterations were performed. Other 
details of the parameters for each simulation are given in 
the figure captions. There are no-slip velocity boundaries 
on the walls. 

The results in Figure 1 aim to compare the cases where 
the director at the wall is pinned or free to rotate. Con- 
sider first the case where the director is free. As a bench 
mark we show in Figure la results from simulations where 
the back-flow, ie the effect of the director field on the flow 
field, is turned off. The fiow profile is as expected the 
parabola of Poiseuille flow. The cholesteric helix drifts 
at a slightly (^ 10%) smaller velocity than the flow. It 
is slightly bent by the fiow into a chevron structTirc. 

The origin of the chevron pattern can be identified by 
considering Eq. 3. Focussing on the center of the chan- 
nel, S is zero, and the term u ■ VQ - typical of this flow 
mode, in which the director field is not constant along the 
fiow direction - must be balanced by a drift of the layers, 
dtQ. However, due to the parabolic shape of the veloc- 
ity field, these cannot cancel exactly and it is necessary 
to allow for a non-zero molecular field, resulting in the 
observed bending. The dominant elastic deformations as- 
sociated with this steady state solution are splay-bend. 
The director field also develops a small component along 
the flow direction. This is caused by the shear forces con- 
tained in the tensor S, which is non-zero away from the 



2 



center of the channel. 
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FIG. 1. Velocity (left) and director fields (right) for (a) 
free boundary conditions, no back-flow; (b) free boundaries 
and back-flow; (c) fixed boundaries, no back-fiow; (d) fixed 
boundaries and back-flow. In (b) and (d) we also show the 
secondary flow (dashed line). This is an averaged velocity be- 
cause the velocity field changes slightly along the flow direc- 
tion y due to the helical arrangement of the director field. The 
system studied corresponds to a cell of thickness L — 2.25 /iin, 
and to a cholesteric liquid crystal with p = 1.6 /Ltm, K ~ 6.3 
pN. The rotational viscosity is 71 = 1 Pois, while the ratio 
between the Leslie viscosities 02 and as is azjoii ~ —0.23 
(flow tumbling). 

Figure lb compares the case when the full equations 



of motion, inchiding back-flow, arc solved. The velocity 
is similar in magnitude, and the parabola flattens. This 
effect becomes more pronounced with an increase in the 
system size (Fig. 2). Our results are consistent with the 
picture described in De Gennes and Prost [1] (although 
these authors consider pinned director boundary condi- 
tions) that the velocity profile in a wide system will be 
flat, decaying to zero near the boundaries in a length 
of order the pitch of the cholesteric helix to satisfy the 
no-slip conditions on the velocity. Also in this case the 
cholesteric helix drifts with the flow. A significant effect 
of the back-flow is that there is much less bending of the 
helix, the pitch remaining constant across the system. 




FIG. 2. Velocity profiles as a function of scaled distance 
along z. The pressure difference is scaled so that the velocity 
field of the liquid crystal in the isotropic phase is the same 
for all three sizes. The other parameters are as in Fig. 1. 

In the steady state we also find a secondary flow, of 
odd parity in z, along the x direction, which attains its 
maxima (in absolute value) close to the boundaries, and 
is zero in the center of the channel. The magnitude of this 
flow is comparable to the maximum fluid velocity along 
y and it gives rise to shear forces which allow Eq. 3 to 
be balanced with a smaller director field deformation. 

Now consider the case where the director is pinned at 
the wall (in such a way that there is no frustration in the 
helix when the system is at rest). Figure Ic shows the 
case without back-flow. The velocity field is, as expected, 
identical to that in Figure la because the director is hav- 
ing no effect on the flow. The helical texture is however 
much more deformed with a substantial component along 
y. This occurs because the director conflguration is un- 
able to drift with the flow and thus the term u ■ VQ (in 
Eq. (3)) requires a larger balancing molecular field. 

When back-flow is turned on for this case there is a 
striking difference in the velocity profile. The net veloc- 
ity is zero, to within the accuracy of the simulations (as 
would be expected from approximate analytic treatments 
[1,4].) Secondary flow is important: it is bigger than the 
maximum velocity attained in the primary flow. Because 
the velocities in the system are very small the director 
field remains very close to its zero-flow configuration. 

We now consider what happens as the pressure differ- 
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cncc is increased. First the chevrons gradually become 
more bent until a threshold forcing at which a new struc- 
ture appears. This structure, shown in Figure 3, has a 
flow-induced twist in the 2;-direction. (The solution pic- 
tured in Figure Ic corresponds to the threshold between 
the chevron and the double twisted texture.) The period 
of the twist along z is roughly equal to the natural twist 
of the cholesteric liquid crystal. This double twist struc- 
ture was found for both flow-aligning (typically found in 
display devices) and flow tumbling (used in some rheo- 
logical experiments [9]) viscosity regimes. 
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Director field for a system with fixed boundary 

Parameters are L = 3 /im, p = 1.6 /im, K ~ 6.3 
~ 0.08 (flow aligning). The maximum velocity in 
equilibrium is ~ 45 pm s~^ (for 71 = 1 Pois). Note the double 
twist texture near the center of the channel. 
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is stabilized by a large enough forcing. 

Our results suggest that in real experiments [5,6,9], for 
slow flows, the helix is likely pinned to its initial config- 
uration (e.g. by irregularities), so that the situation is 
close to that in Figure Id and the viscosity is very large. 
If the pressure difference is increased, either the pinning 
is destroyed or a double twist forms. In both cases the 
viscosity would decrease as observed experimentally. 

In conclusion, we have presented lattice Boltzmann 
simulations able to successfully simulate Poiseuille flow 
in cholesteric liquid crystals in the permeation mode. For 
weak forcing, if the cholesteric helix is pinned at the 
boundaries, we find a remarkable increase in the viscos- 
ity. With free boundaries there is a plug-like velocity 
profile in which the helix drifts with the flow. Beyond a 
threshold, dependent on system parameters, the flow dis- 
torts the helix much more giving rise to a doubly twisted 
director pattern. This approach could also be used to in- 
vestigate, for the flrst time, the rheology of blue phases. 
Predictions can tested with present day microchannel ex- 
periments. In particular a fast Poiseuille flow may pro- 
vide a novel method of inducing double twist in a liquid 
crystal texture thus allowing controlled experiments on 
blue phases in cholesterics. 

This work was supported by EPSRC grant no. 
GR/R83712/01. 



The maximum fluid velocity beyond which double 
twist appears depends sensitively on the system param- 
eters. It is consistently smaller if fixed boundary condi- 
tions are used. For example, for K ^ ^ pN, L ^ 2.25fim, 
a-ija-i ~ —0.23 and p ~ 1.6 m, with fixed boundary 
conditions the velocity threshold is ~ 30/ims~^ while for 
free boundaries it is an order of magnitude bigger. As L 
increases, the threshold velocity decreases if fixed bound- 
aries are used. This effect is somewhat less pronounced 
with free boundaries. The crossover between chevrons 
and double twist is smooth; even at weak forcing one can 
identify a small twist deformation in the z direction. 

It is possible to explain qualitatively the appearance of 
a double twist by again looking at Eq. 3. Let us consider 
the case of anchored boundaries. Near the center of the 
channel, the solution is still determined by balancing the 
molecular field with u- VQ, but since the flow is faster, 
the bending of the chevrons become progressively more 
enhanced. When the deformation is such that two re- 
gions ordered along x are a distance ~ p apart along the 
doubly twisted state is expected to be more sta- 
ble. Such a texture is reminiscent of the doubly twisted 
cylinders which are the basic constituents of blue phases 
[1]. Indeed, in the case of no back-flow, where the simu- 
lations are stable for faster flows, a stack of such doubly 
twisted cylinders, separated by a regular array of defects, 
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